MEMORANDUM  REPORT  NO.  1840 


U.S.  ARMY  BALLISTIC  RESEARCH  LABORATORIES 
ABERDEEN  PROVING  GROUND,  MARYLAND 


6  April  1967 


ERRATA  SHEET 
for 

BRL  Memorandum  Report  No.  1840  entitled  "Optical  Method  for  Analysis  of 
Atmospheric  Effects  on  Laser  Beams",  dated  July  1967* 


Page  16,  line  9.  Change  5*0  to  0.50 

Page  18,  change  the  following  entries  in  Table  II  to  read; 

*  a 

^(calculated) 

Afternoon  0. 38 

2.0  x  10"T  m”1/3 

Night  O.lh 

7.5  x  10"8  m"1/3 

3.  Appendix  -  Re  Figures  A-l  to  A-4,  A-lU:  The  photographs  of  the  beam 
cross  section  should  he  oriented  to  shov  the  step  wedge  along  the 
•left  margin  rather  than  the  bottom  margin.  See  text  (pages  12-13) 
for  optical  system  description. 


BALLISTIC 


RESEARCH  LABORATORIES 


MEMORANDUM  REPORT  NO.  lBkO 
JULY  1967 


OPTICAL  METHOD  FOR  ANALYSIS  OF 
ATMOSPHERIC  EFFECTS  ON  LASER  BEAMS 


Paul  H.  Deitz 

Ballistic  Measurements  Laboratory 


This  material  was  presented  at  the  Symposium  on 
"Modern  Optics,"  sponsored  by  Brooklyn  Polytechnic 
Institute  and  held  at  Waldorf-Astoria  Hotel,  New  York 
City,  21-24  March,  19 67. 


Distribution  0*  this  document  is  unlimited. 


RDT&E  Project  No.  1P523801A286 


ABERDEEN  PROVING  GROUND,  MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 


MEMORANDUM  REPORT  NO.  l840 


PHDeitz/sjv 

Aberdeen  Proving  Ground,  Md. 
July  1967 


OPTICAL  METHOD  FOR  ANALYSIS  OF 
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ABSTRACT 

This  report  describes  the  design  of  an  instrumentation  system  for 
the  study  of  the  effects  of  atmospheric  turbulence  on  a  collimated  laser 
beam  under  near-earth  conditions. 

The  instrumentation  consists  of  a  helium-neon  laser  with  optical 
collimator  and  a  receiving  system  of  24-inch  aperture  with  narrow  band- 
peas  filter. 

The  Appendix  presents  examples  of  beam  cross  section  patterns  for 
different  propagation  conditions.  The  method  of  analyzing  spatial 
intensity  distributions  of  these  patterns  is  described. 
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INTRODUCTION 


Nearly  all  optical  systems  involve  the  transmission  of  electro¬ 
magnetic  radiation  through  an  atmosphere.  For  most  systems  operating  in 
controlled  atmospheres,  the  medium  is  assumed  to  he  homogeneous  and  its 
effect  on  system  performance  negligible.  However,  for  long  paths 
through  the  atmosphere  the  medium  is  not  homogeneous  and  the  index  of 
refraction  varies  along  the  beam;  this  variation  can  be  characterized 
statistically.  Because  of  the  effect  of  the  medium  on  the  transmitted 
beams,  the  performance  of  many  optical  systems  is  finally  limited  by 
atmospheric  conditions. 

To  study  these  effects,  the  Ballistic  Research  Laboratories  (BRL) 
have  designed  and  built  a  high  resolution,  high  stability  optical 
receiving  system.  It  is  designed  specifically  for  the  study  of  light 
propagation  in  the  near-earth  atmosphere.  A  schematic  drawing  of  the 
system  is  shown  in  Figure  1.  This  report  describes  the  design  of  the 
system  and  the  method  used  to  analyze  atmospheric  effects  on  a  collimated 

light  beam  recorded  by  the  system.  Results  of  the  analysis  show  that  the 

i* 

effects  correlate  well  with  those  predicted  hy  Tatarski.  This  system 
should  be  a  valuable  aid  to  a  better  understanding  of  atmospheric  effects 
on  optical  systems . 


Figure  1.  Ray  diagram  of  optical  receiver 

n 

Superscript  markers  denote  references  which,  may  be  found  on  page  2C. 
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Thie  light  source  is  a  helium-neon  laser  with  a  power  of  10-milliwatts . 
The  output  can  he  attenuated  hy  a  polaroid  filter  and  monitored  at  the 
".  laser;-  The  attenuated  beam  is  sent  into  a  collimator  with  a  spatial 
:  filter  located  at  the  focus  of  the  eyepiece  to  remove  intensity  variations 
:  across  the  laser  beam.  The  2-inch  objective  lens  is  defocused  slightly 
so  that  the  transmitted  Airy  disc  is  somewhat  larger  than  the  receiver 
aperture  at  the  range  of  5*+9  meters.  The  laser  and  related  optics  are 
r  held  on  an  aimable  mount.  This  system  is  located  in  a  small  trailer, 
solidly  ground  supported. 

^The  optical  system  is  contained  in  a  25-foot  trailer  making  possible 
easy  transportation  from  one  site  to  another.  On  location,  the  heavily 
reinforced  van  is  jacked  free  from  its  wheels.  Inside  the  trailer,  the 
optical  system  is  mounted  on  an  8-inch  I-beam  frame.  This  frame  is 
lifted  free  of  the  trailer  floor  by  means  of  a  three-point  jacking  system 
supported  directly  on  the  ground. 

Collimated  light  impinging  on  a  parabolic  mirror  2  feet  in  diameter 
with  a  10-foot  focal  length  is  reflected  off-axis  toward  a  plane  mirror 
and  then  a  focus.  Ac  aperture  is  placed  at  this  focus  which  restricts 
the  acceptance  angle  of  the  system  to  about  two  tenths  of  a  degree. 

Vest,  an  f/4,  6-inch  focal  length  lens  is  placed  one  focal  length  from 
the  aperture.  Thus,  the  emerging  light  is  collimated,  making  possible 
Introduction  of  an  interference  filter  which  eliminates  nearly  all  of  any 
remaining  background  illumination.  Beyond  this  first  lens  and  filter  is 
an  image  of  the  spatial  intensity  distribution  at  the  2-foot  parabolic 
mirror.  Alongside  this  image,  an  illuminated  step  wedge  is  projected  in 
conjunction  with  fiducial  lights  to  enable  frame -to-frame  correlation 
between  a  timed  series  of  photographs.  Both  the  intensity  distributions 
of  the  parabolic  mirror  and  the  step  wedge  then  are  imaged  by  a  second 
lens  to  a  sensor  plane  at  a  magnification  compatible  with  a  chosen 
format. 
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CALIBRATION 

An  electromagnetic  wave  is  completely  characterized  by  its  amplitude 
end  phase.  The  optical  receiver  was  set  up  to  record  the  amplitude  of 
the  impinging  radiation  in  terms  of  its  related  intensity.  A  35-milIi- 
meter  motion  picture  camera  was  placed  in  the  film  plane  of  the  receiver 
to  record  the  spatial  intensity  pattern  of  the  incident  radiation  along 
with  the  calibrated  step  wedge.  The  wedge  was  calibrated  at  the  film 
plane  to  an  accuracy  of  5  percent.  The  intensity  variation  across  the 
extent  of  the  entrance  aperture  of  the  receiver  was  measured  by  illumi¬ 
nating  the  parabolic  mirror  with  the  far-field  pattern  and  scanning  the 
image  plane.  Over  the  central  18  inches  of  the  mirror,  the  intensity 
variation  was  less  than  6  percent.  A  geometrical  distortion  and  resolu¬ 
tion  test  ’.as  made  by  placing  a  test  screen  in  front  of  the  receiver  and 
back-illuminating  it  with  the  divergent  laser  beam.  A  grid  array  of 
holes  on  10-centimeter  centers  was  drilled  in  the  board  of  «mch  size  as 
to  give  a  2  millimeter  diameter  Fraunhofer  diffraction  pe.t\  ■  on  the 
parabolic  mirror  surface.  Pictures  were  taken  ~>t  this  pattern  and 
measured  to  a  1-micron  precision  on  a  comparator.  Accurate  measurements 
were  made  of  the  grid  array,  and  a  least  squares  fit  computer  program  was 
run  on  these  data.  The  standard  deviation  of  the  porition  of  the  image 
points  related  to  object  space  and  compared  with  the  grid  array  was  0.7 
millimeter  with  a  resolution  of  0.25  millimeter  or  better  for  the  total 
receiver.  A  recently  modified  fiducial  mark  projection  system  gives 
positional  accuracy  in  the  image  plane  to  within  10  microns.  In  object 
space,  this  corresponds  to  a  positional  accuracy  of  1.3  millimeters. 


ANALYTIC  METHOD 


The  type  of  optical  data  being  sought  initially  was  beam  intensity 
as  a  function  of  position  on  a  horizontal  and  a  vertical  diameter  of  the 
parabolic  mirror.  Various  films  had  been  tested  to  determine  the 
dynamic  range  and  gamma  of  the  linear  portion  of  the  characteristic  curve 
of  the  film.  For  each  type  of  film  and  exposure  time,  a  particular 
absolute  intensity  variation  of  the  step  wedge  covers  the  dynamic  range 
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or  usable  portion  of  the  film  curve.  One  of  the  brighter  steps  gives  an 
intensity  at  the  top  of  the  dynamic  range  of  the  film.  Heutral  density 
filters  were  inserted  at  the  optical  receiver  until  the  most  intense 
radiation  corresponds  to  the  intensity  of  the  brightest  usable  step. 

After  the  data  were  taken  and  the  film  developed*  the  individual  frames 
were  scanned  consecutively  on  a  densitometer.  Tue  scans  were  always  over 
the  same  2  diameters  of  the  parabolic  mirror.  The  densities  as  a  function 
of  position  were  recorded  and  then  converted  to  relative  intensities 
through  the  characteristic  curve  of  the  film  obtained  from  the  s'tep 
wedge. 

A  series  of  trial  photographs  were  made  to  determine  the  feasibility 
of  the  technique.  The  path  length  is  5*+9  meters,  averaging  about  2  meters 
above  short  grass.  The  location  was  the  Electromagnetic  Propagation  Range 
at  the  BRL.  The  range  features  an  extensive  array  of  meteorological 
sensors  with  an  electronic  data  acquisition  system  (Figure  2). 


Figure  2.  Electromagnetic  propagation  range  with 
30-foot  sensor  stations 

Figure  A-l  in  the  Appendix  shows  a  mottled  pattern  typical  for  a 
day  with  high  scintillation.  The  step  wedge  extends  vertically  in  the 
optical  system.  Since  the  system  incorporates  an  optical  flat,  there 
is  a  horizontal,  but  not  a  vertical,  image  inversion.  Thus,  the  image 
can  be  understood  to  he  the  intensity  distribution  that  would  be  seen  if 
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a  2-foot  section  of  ground  glass  were  to  be  substituted  for  the  parabolic 
mirror  with  the  observer  located  behind  the  ground  glass  looking  back 
toward  the  source.  This  photograph  was  taken  in  the  susner  during  a 
mid-afternoon.  The  exposure  time  for  all  of  the  photographs  is  1/500 
second.  There  was,  therefore,  no  image  blur  in  the  horizontal  direction. 
Figure  A-2  is  a  photograph  taken  at  night  under  medium  scintillation 
conditions.  Since  the  atmosphere  was  fairly  quiet,  the  high  degree  of 
coherence  of  the  laser  was  maintained,  resulting  in  the  characteristic 
diffraction  effects  from  dust  particles  on  the  lenses. 

Figure  A- 3  shows  an  example  of  a  fleeting  optical  effect  observed 
after  dusk  one  summer  evening.  The  pattern  consisted  of  large,  slowly 
changing  cells  which  built  up  gradually,  lasting  in  this  condition  for 
about  a  minute.  The  pattern  broke  up  quickly  into  much  higher  spatial 
frequencies;  once  more  built  up  to  the  large  cell  structure  for  a 
minute;  and  finally  broke  up  altogether. 

Figure  A-U  shows  one  of  two  photographs  of  which  a  quantitative 
analysis  has  been  made.  The  step  wedge  at  the  left  of  the  photograph  was 
scanned  and  a  computer-fi ited  characteristic  curve  derived  (Figure  A-5). 
The  relative  densities  are  plotted  vertically  and  the  logarithm  of  the 
exposure  is  shown  to  the  right.  The  antilogs  are  also  given.  The 
extremes  of  the  exposure  are  indicated  for  the  horizontal  scan  through 
the  mirror.  Using  this  curve,  the  output  of  the  densitometer  was 
computer-corrected  to  relative  intensities,  and  plotted  (Figure  A-6) . 

THEORY  AND  RESULTS 

Tatarski1  predicts  the  two-dimensional  spatial  spectral  densities 
in  a  plane  perpendicular  to  the  direction  of  propagation  for  a  plane, 
monochromatic  wave.  The  spatial  spectral  densities  are  postulated  to  be 
isotropic  and  to  be  a  function  tally  of  the  wavelength,  X,  and  the  path- 
length,  L.  For  the  case  where  the  wind  direction  angle  relative  to  the 
optical  path  is  much  larger  than  the  quantity  (A/L) 9  the  intensity 
structure  is  assumed  to  be  ”frozen-inn  and  carried  across  the  optical 
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path  by  the  normal  wind  component ,  vfi.  Tatarski  thus  translates  from  a 
two-dime*  sional  spatial*  spectral  density  domain  to  a  one-dimensional 
vc^tdjie-  domain .  •' 

Tatarskiis  experiment,  a  weakly  divergent  incandescent  beam  is 

propagated  over  a  homogeneous  path  within  2  meters  of  the  ground.  The 

received  light  is  monitored  by  a  photocell  with  an  aperture  smal  1  »nmigh 

,/  (2  millimeters)  to  resolve  the  hipest  spatial  frequencies.  As  the  normal 

wind  component  takes  the  turbulent  medium  across  the  beam  path,  a  one- 

'  dimensional  intensity  profile  is  recorded.  Since  the  optical  path  is 

near  the  ground,  the  vertical  wind  component  is  necessarily  small,  so 

that  the  profile  is  predominantly  horizontal.  The  magnitude  of  vq  is 

monitored  during  the  experiment.  Figure  A-7  shows  the  results  of  such 

measurement.  The  output  of  the  photocell  is  sent  into  a  time  spectrum 

analyzer.  The  power,  W(f),  is  multiplied  by  temporal  frequency  and 

plotted  versus  the  log^  of  the  temporal  frequency.  As  can  be  expected, 

as  v^  goes  to  higher  values,  the  peak  power  point  of  the  curve  also  moves 

.to  higher  frequencies.  Next,  all  these  curves,  along  with  data  taken  at 

different  pathlengths,  are  normalized  to  a  v  of  1  meter  per  second  and 

n 

a  path  of  1,000  meters  (Figure  A-8) .  The  curves  then  assume,  generally, 
the  same  shape. 

Tatarski *8  time  spectrum  is  related  to  the  space  spectrum  by 

T  »  S  •  vn  (1) 

where  the  space  spectrum,  S,  is  measured  in  cycles  per  meter.  When  the 
space  spectrum  is  thus  measured,  it  is  equivalent  to  the  time  spectrum 
inherently  normalized  to  e.  vq  of  1  meter  per  second.  To  test  this 
relationship,  a  power-density,  spectrum  analysis  vas  made  using  the 
horizontal  intensity  profile  shown  in  Figure  A- 6.  Figure  A-9,  then,  shows 
the  relative  powers  as  a  function  of  spatial  frequencies.  Hext,  the  power 
was  multiplied  by  frequency  and  plotted  on  a  logarithmic  abscissa 
a  (Figure  A-1C) .  Since  the  spectra  have  frequency  plotted  in  cycles  per 
meter  and  are  thus  normalized  for  a  normal  wind  component  of  1  meter  per 
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second,  the  space  spectra  units  convert  directly  to  time  spectra  units. 
The  roughness  of  the  curve  is  indicative  of  the  shortness  of  the  sample 
length.  To  get  a  more  definitive  plot,  many  such  samples  should  be 
analyzed  and  averaged.  In  order  to  approximate  the  effect  of  a  longer 
sample  length,  a  smoothing  program  was  applied  to  Figure  A-10  resulting 
in  the  curve  shown  in  Figure  A-ll.  This  reveals  a  much  greater  likeness 
to  the  Tatarski  time  plots  but  the  spatial  spectra  were  normalized  by 
different  schemes  and  the  ordinates  cannot  be  compared.  With  this  method 
of  data  acquisition  there  is  no  time  lag  from  the  beginning  of  the 
intensity  scan  to  the  end  so  no  assumption  about  the  degree  to  which  the 
turbulent  cell  structure  is  ”frozen-in"  is  necessary  to  infer  the  nature 
of  the  true  spatial  frequency  characteristics .  Under  the  assumption  that 
the  photograph  was  indeed  isotropic,  zhe  unsmoothed  vertical  and  hori¬ 
zontal  power  density  spectra  were  averaged,  multiplied  by  frequency,  and 
plotted  on  a  logarithmic  abscissa  (Figure  A-12).  These  combined  spectra 
units  were  then  smoothed  (Figure  A-13),  and  it  was  this  curve  that  was 
used  for  comparison  to  the  theory. 

By  the  same  process,  a  photograph  taken  during  a  night  propagation 
run  was  also  analyzed  (Figure  A-lk).  As  in  the  earlier  photograph  taken 
at  night.  Figure  A-2,  the  effect  of  highly  coherent  light  operating  on 
the  optical  receiver  can  be  teen  as  diffraction  rings  in  the  image  plane. 
Figure  A-15  shows  the  power  density  spectrum,  and  Figures  A-l6  to  A-19 
show  the  frequency  plots  smoothed  and  averaged  as  before. 

Tatarski  predicts  that  the  frequency  at  which  the  maximum  power 
occurs  is  given  by  the  formula 


(2) 


where  f  is  defined  as  the  average  of  the  two  frequencies  occurring  at 
w 

one-half  the  peak  power.  Theoretically,  A  is  given  as  0.55,  while 
experimentally  Tatarski  obtains  a  value  of  0.3 2.  Using  Figures  A-13  and 
A-19,  an  average  A  of  0.222  was  obtained  for  the  two  photographs. 
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Although  not  explicitly  stated  later  in  Tatar ski,  f  is  derived  for  a 

m 

plane  wave  source.  The  difference  between  the  theoretical  and  experi¬ 
mental  values  is  not  explained.  If,  in  a  very  simple  model,  all  of 
the  diffraction  is  assumed  to  take  place  at  a  plane  located  midway  between 
the  transmitter  and  receiver,  then  the  spatial  frequencies  will  be  scaled 
down  at  the  receiver  by  the  ratio  of  the  beam  width  at  the  path  midpoint 
to  the  beam  width  at  the  path  terminus.  As  transmitted  light  is  changed 
from  a  collimated  beam  to  a  point  source,  B  goes  from  a  value  of  1.0  to 
5.0.  Table  I  shows  this  ’'divergence  factor"  for  both  experiments.  This 
then  is  to  be  applied  as  a  correction  factor  to  the  coefficient  0.55  for 
a  plane  wave.  The  product  of  these  two  is  in  the  column  titled  "A, 
calculated."  This  value  can  be  compared  with  the  measured  value  of  A. 


Table  I 


.  _  „  .  width  of  beam  at  midooint  ..  _  . 

A  -  B(0.55)  B  -  riath  of  tern  st  temlmia  *  ^ 


B 


calculated 


calculated 


A 


measured 


Plane  Wave 

1.0 

0.55 

— 

Tatarski 

0.6 

0.33 

0.32 

BRL 

0.50h 

0.277 

0.222 

The  importance  of  resolving  the  highest  spatial  frequencies  in 

object  space  was  inferred  earlier.  The  influence  of  the  diameter  of  the 

photocell  aperture  on  the  frequency  distribution  has  been  predicted  by 

2 

Tatarski  and  measured  by  others.  With  the  present  demagnification  from 
object  space  to  the  image  plane  and  the  size  of  the  scanning  aperture  of 
the  densitometer,  the  effective  aperture  size  in  object  space  is  1.26 
millimeters.  To  show  the  effect  of  an  increasing  one-dimensional  aper¬ 
ture,  the  intensities  shown  in  Figure  A-6  were  averaged  over  windows 
with  widths  of  5>  9,  17 »  and  33  samples.  Each  of  the  smoothed  scans  was 
subjected  to  frequency  analysis,  smoothed,  and  plotted  (Figure  A-20). 

The  smoothed  frequency  times  power  plot  of  Figure  A-ll  is  also  shown  for 


16 


comparison.  The  effective  sizes  of  the  apertures  in  increasing  widths 
are  1.26,  5*16,  9.07,  16.9,  and  32.5  millimeters .  As  the  widths  are 
increased,  both  the  decrease  and  shift  in  peak  power  and  the  loss  of 
high  frequency  response  are  readily  apparent. 

The  shape  and  maxima  of  the  spatial  frequency  spectra  are  not 
supposed  to  vary  for  a  given  pathlength  or  wavelength.  However,  the 
intensities,  which  are  postulated  to  he  log  normally  distributed,  are 
predicted  to  exhibit  a  variance  proportional  to  the  fluctuations  in  the 
index  field.  Tatarski's  formula  giving  this  relationship  is 


o2  -  1.23  0  2  k7/6  Ln/6 

n 


(3) 


where  a  is  the  standard  deviation  of  the  log  normal  distribution,  CQ  is 
the  index  structure  constant,  k  is  the  wave  number  of  the  light,  and  L 
is  the  pathlength.  To  test  this  relationship,  the  horizontal  intensities 
from  Figures  A-h  and  A~lU  were  each  normalized  by  division  of  their 
average  intensity.  The  percent  number  of  occurrences  was  then  plotted 
versus  the  normalized  intensities  to  form  distribution  curves  (Figures 
A-21  and  A-22).  Figure  A-21  from  the  afternoon  propagation  data  shows 
a  much  greater  skewness,  indicating  higher  scintillation  conditions. 

Next,  the  integral  of  each  of  these  curves  was  taken  to  yield  the 
cumulative  distribution  curves.  The  integrated  percent  occurrences  were 
plotted  on  probability  paper  versus  the  log  of  the  normalized  intensities 
(Figures  A-23  and  A-2h).  Both  sets  of  data  conform  closely  to  log  normal 
distributions.  The  divergence  at  the  extremes  of  the  curve  in  Figure  A-23 
can  be  interpreted  in  the  following  way.  The  first  point  was  taken  from 
the  very  edge  of  the  mirror  where  the  calibration  has  indicated  a  sharp 
attenuation  of  intensities  due  to  vignetting.  Ifce  intensities  at  the 
upper  end  of  the  curve  came  from  the  top  of  the  characteristic  curve. 

Since  they  were  forcing  the  dynamic  range  of  the  film,  their  true  values 
were  likely  higher  than  those  indicated.  Although  no  CQ  was  measured 
for  the  two  caseB,  the  cumulative  distribution  curves  were  used  to  derive 
the  standard  deviation.  Equation  (3)  vas  solved  for  Cn,  and  the  appropri¬ 
ate  values  were  inserted;  Table  II  lists  the  results. 
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Table  II 


0 

C  (calculated) 
n 

Afternoon 

0.852 

1.96  x  10'T  nf1/3 

Night 

0.332 

7.63  x  10"8  m“1/3 

^Weak  turbulence 

C  =  8  x  10*9  m-1^3 
n 

Intermediate  turbulence 

C  *  k  x  10-8  a"1^3 
n 

Strong  turbulence 

Cq  =  5  x  10“7  m"1/3 

•3 

Coatparing  the  calculated  values  of  with  those  given  bj  Divis  for 
three  representative  cases  indicates  reasonable  agreement  for  the  inter¬ 
mediate  to  strong  scintillation  conditions. 

SUMMARY  AND  CONCLUSION 

Initial  data  analysis  of  the  spatial  intensity  distributions  from 
the  optical  receiver  shows  a  high  degree  of  correlation  to  the  theoreti¬ 
cal  predictions  of  Tatarski.  When  the  divergence  of  the  transmitted  beam 
is  taken  into  consideration,  the  predicted  frequency  of  maximum  power 
compares  favorably  with  the  measured  value.  Widening  of  the  sampling 
aperture  is  shown  both  to  attenuate  the  high  frequency  resolution  and  to 
shift  the  peak  power  point  to  lower  frequencies.  Calculated  values  of 
the  index-structure  constant  from  the  measured  variance  of  the  intensity 
profiles  conform  to  the  proper  orders  of  scintillation  magnitude. 

During  normal  data  runs,  a  time  series  of  photographs  is  taken  to 
make  the  results  statistically  meaningful.  Since  in  this  particular 
approach  the  scanning  is  not  a  function  of  the  normal  wind  component, 
any  intensity  profile  can  be  examined,  including  the  vertical,  allowing 
a  check  of  the  assumption  of  optical  isotropy  both  by  electronic  reduction 
and  direct  optical  correlation  techniques.  By  this  method,  the  spatial 
frequency  domain  may  be  converted  to  a  temporal  frequency  dome,  .n  which 
is  inherently  normalized  for  a  wind  component  normal  to  the  direction  of 


beaa  propagation  of  1  meter  per  second.  No  assumption  about  the  degree 
to  which  the  turbulent  cell  structure  is  "frozen-in"  is  necessary  to 
infer  the  nature  of  the  spatial  frequency  statistics.  In  addition,  the 
high  resolution  and  good  signal-to-noise  characteristics  of  the  receiver 
make  accurate  measurements  of  beam  pointing  reliability  possible  for 
various  system  applications. 

This  BRL  approach  will  help  to  give  a  better  quantitative  under¬ 
standing  of  the  constraints  placed  on  nearly  all  optical  systems  working 
in  the  atmosphere. 
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Figure  A-1.  Beam  cross-section  photograph  taken  under  high 

scintillation  conditions 
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Figure  A-2.  Beam  cross-section  photograph  taken  under  medium 

scintillation  conditions 


Figure  A-4.  Sample  of  beam  cross-section  data  used  for  quantitative 
analysis  (afternoon  propagation) 
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Figure  A-7,  Temporal  frequency  times  spectral  density  vs  temporal  frequency 
(Tatarskl,  Wave  Propagation  in  a  Turbulent  Medium.  McGraw-Hill,  1961,  page  220) 


Figure  A-8.  Temporal  frequency  times  spectral  density  vs  temporal  frequency 
normalized  to  a  normal  wind  component  of  1  meter/second 

(Tartarskl,  Wave  Propagation  In  a  Turbulent  Medium,  McGraw-Hill,  1961,  page  222) 
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Figure  A-9. 
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•  power  spectral  density  plot  of  horizontal  'intensity 
profile  (afternoon  propagation) 
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Figure  A-13.  Spatial  frequency  V  as  averaged  smoothed  horizontal  and 
vertical  power  density  spectra  vs  spatial  frequency 
(afternoon  propagation) 
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Figure  A-14.  Sample  of  beam  cross-section  data  used  for  quantitative 

analysis  (night  propagation) 
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Figure  A- 18.  Spatial  frequency  times  averaged  horizontal  and  vertical 
power  density  spectra  vs  spatial  frequency 
(night  propagation) 
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gure  A-20.  Spatial 
spatial  freque 
apertures 
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